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Abstract 


We present an experiment in which a one-bit memory is con¬ 
structed, using a system of a single colloidal particle trapped in a 
modulated double-well potential. We measure the amount of heat dis¬ 
sipated to erase a bit and we establish that in the limit of long erasure 
cycles the mean dissipated heat saturates at the Landauer bound, i.e. 
the minimal quantity of heat necessarily produced to delete a classical 
bit of information. This result demonstrates the intimate link between 
information theory and thermodynamics. To stress this connection we 
also show that a detailed Jarzynski equality is verified, retrieving the 
Landauer’s bound independently of the work done on the system. The 
experimental details are presented and the experimental errors care¬ 
fully discussed 
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1 A link between information theory and 
thermodynamics 

The Landauer’s principle was first introduced by Rolf Landauer in 1961 [10]. 
It states that any logically irreversible transformation of classical information 
is necessarily accompanied by the dissipation of at least kpT In 2 of heat per 
lost bit, where fee is the Boltzmann constant and T is the temperature. This 
quantity represents only ~ 3 x 10 -21 J at room temperature (300 K) but is 
a general lower bound, independent of the specific kind of memory system 
used. 

An operation is said to be logically irreversible if its input cannot be 
uniquely determined from its output. Any Boolean function that maps sev¬ 
eral input states onto the same output state, such as AND, NAND, OR and 
XOR, is therefore logically irreversible. In particular, the erasure of informa¬ 
tion, the RESET TO ZERO operation, is logically irreversible and leads to 
an entropy increase of at least k b In 2 per erased bit. 

A simple example can be done with a 1-bit memory system (i.e. a systems 
with two states, called 0 and 1) modelled by a physical double well potential 
in contact with a single heat bath. In the initial state, the bistable potential 
is considered to be at equilibrium with the heat bath, and each state (0 
and 1) have same probability to occur. Thus, the entropy of the system is 
S = kp In 2, because there are two states with probability 1/2. If a RESET 
TO ZERO operation is applied, the system is forced into state 0. Hence, 
there is only one accessible state with probability 1, and the entropy vanishes 


2 







5 = 0. Since the Second Law of Thermodynamics states that the entropy 
of a closed system cannot decrease on average, the entropy of the heat bath 
must increase of at least k B bi 2 to compensate the memory system’s loss of 
entropy. This increase of entropy can only be done by an heating effect: the 
system must release in the heat bath at least k B T In 2 of heat per bit erased 1 . 

For a reset operation with efficiency smaller than 1 (he. if the operation 
only erase the information with a probability p < 1), the Landauer’s bound 
is generalised: 


( Q) > k B T [In 2 + pln(p) + (1 -p)ln(l - p)\ (1) 

The Landauer’s principle was widely discussed as it could solve the para¬ 
dox of Maxwell’s “demon” [2,4,14], The demon is an intelligent creature 
able to monitor individual molecules of a gas contained in two neighbouring 
chambers initially at the same temperature. Some of the molecules will be 
going faster than average and some will be going slower. By opening and 
closing a molecular-sized trap door in the partitioning wall, the demon col¬ 
lects the faster (hot) molecules in one of the chambers and the slower (cold) 
ones in the other. The temperature difference thus created can be used to run 
a heat engine, and produce useful work. By converting information (about 
the position and velocity of each particle) into energy, the demon is therefore 
able to decrease the entropy of the system without performing any work him¬ 
self, in apparent violation of the Second Law of Thermodynamics. A simpler 
version with a single particle, called Szilard Engine [20] has recently been 
realised experimentally [21], showing that information can indeed be used to 
extract work from a single heat bath. The paradox can be resolved by noting 
that during a full thermodynamic cycle, the memory of the demon, which is 
used to record the coordinates of each molecule, has to be reset to its initial 
state. Thus, the energy cost to manipulate the demon’s memory compensate 
the energy gain done by sorting the gas molecules, and the Second Law of 
Thermodynamics is not violated any more. 

More information can be found in the two books [11,12], and in the very 
recent review [13] about thermodynamics of information based on stochastic 
thermodynamics and fluctuation theorems. 

In this article, we describe an experimental realization of the Landauer’s in¬ 
formation erasure procedure, using a Brownian particle trapped with optical 
tweezers in a time-dependent double well potential. This kind of system was 
theoretically [19] and numerically [6] proved to show the Landauer’s bound 

1 It it sometimes stated that the cost is k B T hi2 per bit written. It is actually the same 
operation as the RESET TO ZERO can also be seen to store one given state (here state 
0), starting with an unknown state. 
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Figure 1: Schematic representation of the one-bit memory system, made of 
one particle trapped in a double well potential. 


k-^T \a.2 for the mean dissipated heat when an information erasure opera¬ 
tion is applied. The results described in this article were partially presented 
in two previous articles [3,5], and were later confirmed by two independent 
experimental works [8,15]. 

The article is organized as follows. In section 2 we describe the experi¬ 
mental set-up and the one bit memory. In section 3 the experimental results 
on the Landauer’s bound and the dissipated heat are presented. In section 
4 we analyze the experimental results within the context of the Jarzinsky 
equality. Finally we conclude in section 5. 

2 Experimental set-up 

2.1 The one-bit memory system 

The one-bit memory system, is made of a double well potential where one 
particle is trapped by optical tweezers. If the particle is in the left-well the 
system is in the state “0”, if the particle is in the right-well the system in 
the state “1” ( see fig. 1). 

The particles are silica beads (radius R = 1.00 ± 0.05 pm), diluted a low 
concentration in bidistilled water. The solution is contained in a disk-shape 
cell. The center of the cell has a smaller depth (~ 80 pm) compared to the 
rest of the cell (~ 1 mm), see figure 2. This central area contains less particles 
than the rest of the cell and provides us a clean region where one particle 
can be trapped for a long time without interacting with other particles. 

The double well potential is created using an Acousto-Optic Deflector 
which allows us to switch very rapidly (at a rate of 10 kHz) a laser beam 
(wavelength A = 1024 nm) between two positions (separated by a fixed dis- 
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Figure 2: Schematic representation of the cell used to trap particles dispersed 
in water (view from the side). The central part has a smaller gap than the 
rest of the cell. 


tance d ~ 1pm). These two positions become for the particle the two wells 
of the double well potential. The intensity of the laser / can be controlled 
from 10 mW to more than 100 mW, which enables us to change the height of 
the double well potential’s central barrier 2 . A NanoMax closed-loop piezo¬ 
electric stage from Thorlabs® with high resolution (5 nm) can move the cell 
with regard to the position of the laser. Thus it allows us to create a fluid 
flow around the trapped particle. The position of the bead is tracked using 
a fast camera with a resolution of 108 nm per pixel, which after treatment 
gives the position with a precision greater than 5 nm. The trajectories of the 
bead are sampled at 502 Hz. See figure 3. 

The beads are trapped at a distance h = 25 pm from the bottom of the cell. 
The double well potential must be tuned for each particle, in order to be as 
symmetrical as possible and to have the desired central barrier. The tuning 
is done by adjusting the distance between the two traps and the time that 
the laser spend on each trap. The asymmetry can be reduced to ~ 0.1 k^T. 
The double well potential U 0 (x, /) (with x the position and / the intensity of 
the laser) can simply be measured by computing the equilibrium distribution 
of the position for one particle in the potential: 



( 2 ) 


One typical double well potential is shown in figure 4. 

2.2 The information erasure procedure 

We perform the erasure procedure as a logically irreversible operation. This 
procedure brings the system initially in one unknown state (0 or 1 with same 

2 The values are the power measured on the beam before the microscope objective, so 
the “real” power at the focal point should be smaller, due to the loss in the objective. 
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Figure 3: Schematic representation of optical tweezers set-up used to trap 
one particle in a double well potential. The Acousto-Optic Deflector (AOD) 
is used to switch rapidly the trap between two positions. The NanoMax piezo 
stage can move the cell with regard to the laser, which creates a flow around 
the trapped particle. “M” are mirrors and “DM” is a dichroic mirror. 



Figure 4: Double well potential measured by computing the equilibrium dis¬ 
tribution of one particle’s positions, with I\ aser = 15 mW. Here the distribu¬ 
tion is computed on 1.5e6 points sampled at 502 Hz(ie.a50 minlong measure¬ 
ment). The double well potential is well fitted by a 6 th degree polynomium. 
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probability) in one chosen state (we choose 0 here). It is done experimentally 
in the following way (and summarised in figure 5a): 

• At the beginning the bead must be trapped in one well-defined state 
(0 or 1). For this reason, we start with a high laser intensity (/high = 
48 mW) so that the central barrier is more than 8 k^T. In this situation, 
the characteristic jumping time (Kramers Time) is about 3000 s, which 
is long compared to the time of the experiment, and the equivalent 
stiffness of each well is about 1.5pN/pm. The system is left 4 s with 
high laser intensity so that the bead is at equilibrium in the well where 
it is trapped 3 . The potential £/ 0 (x,/high) is represented in figure 5b 1. 

• The laser intensity is first lowered (in a time 7] ow = 1 s) to a low value 
(/low = 15 mW) so that the barrier is about 2.2 k^T. In this situation 
the jumping time falls to ~ 10 s, and the equivalent stiffness of each well 
is about 0.3pN/pm. The potential Uo(x,I\ ow ) is represented figure 5b 

2 . 

• A viscous drag force linear in time is induced by displacing the cell with 
respect to the laser using the piezoelectric stage. The force is given by 
f = 'yv where 7 = 67 vRrj (?/ is the viscosity of water) and v the speed 
of displacement. It tilts the double well potential so that the bead 
ends always in the same well (state 0 here) independently of where it 
started. See figures 5b 3 to 5. 

• At the end, the force is stopped and the central barrier is raised again 
to its maximal value (in a time T high = Is). See figure 5b 6 . 

The total duration of the erasure procedure is /i ow +T+ Thigh- Since we kept 
Tiow = Thigh = Is, a procedure is fully characterised by the duration r and 
the maximum value of the force applied / max . Its efficiency is characterized 
by the “proportion of success” Ps, which is the proportion of trajectories 
where the bead ends in the chosen well (state 0 ), independently of where it 
started. 

Note that for the theoretical procedure, the system must be prepared in 
an equilibrium state with same probability to be in state 1 than in state 0 . 
However, it is more convenient experimentally to have a procedure always 
starting in the same position. Therefore we separate the procedure in two 
sub-procedures: one where the bead starts in state 1 and is erased in state 0 , 
and one where the bead starts in state 0 and is erased in state 0. The fact that 

3 The characteristic time for the particle trapped in one well when the barrier is high 
is 0.08 s 
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Figure 5: Schematical representation of the erasure procedure. The potential 
felt by the trapped particle is represented at different stages of the procedure 
(1 to 6). For 1 and 2 the potential Uq(x,I) is measured. For 3 to 5 the 
potential is constructed from [/ 0 (x, ii ow ) knowing the value of the applied 
drag force. 


the position of the bead at the beginning of each procedure is actually known 
is not a problem because this knowledge is not used by the erasure procedure. 
The important points are that there is as many procedures starting in state 
0 than in state 1, and that the procedure is always the same regardless of 
the initial position of the bead. Examples of trajectories for the two sub¬ 
procedures 1 —> 0 and 0 —> 0 are shown in figure 6. 
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(a) 1 —> 0 




(b) 0 —>■ 0 


Figure 6: Examples of trajectories for the erasure procedure, t = 0 corre¬ 
sponds to the time where the barrier starts to be lowered. The two possibil¬ 
ities of initial state are shown. 






































3 Landauer’s bound for dissipated heat 


3.1 Computing the dissipated heat 

The system can be described by an over-damped Langevin equation: 

f)TT 

ix = +f{t) + i{t) (3) 

with x the position of the particle 4 , x — its velocity, 7 = 6nRi] the fric¬ 
tion coefficient (77 is the viscosity of water), Uq(x, /) the double well potential 
created by the optical tweezers, f(t) the external drag force exerted by dis¬ 
placing the cell, and £(t) the thermal noise which verifies (£(£)) = 0 and 
= ^k^TSit — t'), where (.) stands for the ensemble average. 
Following the formalism of the stochastic energetics [17], the heat dissi¬ 
pated by the system into the heat bath along the trajectory x(t) between 
time t — 0 and t is: 

Qo,t = [ - (f (0 - 7 Ht')) x(t') d t'. (4) 

Jo 

Using equation 3, we get: 

Qo,t = J (~~^~(x,I) + x(t') dt'. (5) 

For the erasure procedure described in 2.2 the dissipated heat can be decom¬ 
posed in three terms: 


Qe 


Q ba 


^potential “h Cd: 


rag 


( 6 ) 


Where: 


• Cbarrier is the heat dissipated when the central barrier is lowered and 
risen (/ = 0 during these stages of the procedure): 


Q ba 


dUo, 

dx 


(x, I) ) xdt'+ 


f*^iow+ r +Thigh 


Ti ow +r 


dUo 

dx 


(. x,I ) ) xdt' 

(7) 


• ^potential is the heat dissipated due to the force of the potential, during 
the time r where the external drag force is applied (the laser intensity 
is constant during this stage of the procedure): 

/ Tlow'fT" / \ 

y--^-(x,I)J xdt' (8) 

4 We take the reference x = 0 as the middle position between the two traps. 
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• Qdrag is the heat dissipated due the external drag force applied dur¬ 
ing the time r (the laser intensity is constant during this stage of the 
procedure): 

r Tiow+r 


Q d 


rag 


fx dt' 


(9) 


The lowering and rising of the barrier are done in a time much longer 
than the relaxation time of the particle in the trap (~ 0 . 1 s), so they can 
be considered as a quasi-static cyclic process, and do not contribute to the 
dissipated heat in average. The complete calculation can also be done if we 
assume that the particle do not jump out of the well where it is during the 
change of barrier height. In this case, we can do a quadratic approximation: 
Uo(x, I) = —k(I)x where k is the stiffness of the trap, which evolves in time 
because it depends linearly on the intensity of the laser. Then: 

j "low \ r'l low 

(Qlowering) = (J ~kxX dt' ) = J ~~ d ((x 2 )) . (10) 


Using the equipartition theorem (which is possible because the change of 
stiffness is assumed to be quasi-static), we get: 

/•Tw 


(Q lowering) 


, k^T f 1 \ ksT 
-k —-— d ( — | = —— in 
k 


h 


low 


high 


(ii) 


j) = ^ln and it follows directly 


The same calculation gives (Qnsmg/ — 2 

that (d) barrier) (^lowering) T (Qrising) 0. 

The heat dissipated due to the potential when the force is applied is also 
zero in average. Indeed, the intensity is constant and Uq becomes a function 
depending only on x. It follows that: 


(^potential) 


rhow+T 


'how 


dUo 

dx 


dx = 


U 0 (x) 


"I ^(how) 
®(how+T) 


( 12 ) 


Since the potential is symmetrical, there is no change in Uq when the bead 
goes from one state to another, and {U 0 (x(T\ ow )) — U 0 (x(T\ ow + r))) = 0. 

Finally, since we are interested in the mean dissipated heat, the only rele¬ 
vant term to calculate is the heat dissipated by the external drag force: 


fhow+T 


Q drag 


7 v(t')xdt'. 


(13) 


Where 7 is the known friction coefficient, v(t) is the imposed displacement of 
the cell (which is not a fluctuating quantity) and x can be estimated simply: 


x(t + St/2) 


x{t + St) — x(t) 

St 


(14) 
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We measured Qdrag for several erasure procedures with different parameters 
r and / max . For each set of parameters, we repeated the procedure a few 
hundred times in order to compute the average dissipated heat. 

We didn’t measure (^lowering and Qr isi ng because it requires to know the 
exact shape of the potential at any time during lowering and rising of the 
central barrier. The potential could have been measured by computing the 
equilibrium distribution of one particle’s positions for different values of I. 
But these measurements would have been very long since they require to 
be done on times much longer than the Kramers time to give a good es¬ 
timation of the double well potential. Nevertheless, we estimated on nu¬ 
merical simulations with parameters close to our experimental ones that 
(^lowering + ^potential + Qrising) ~ 0.7 k B T which is only 10 % of the Landauer’s 
bound. 

3.2 Results 

We first measured Ps the proportion of success for different set of r and / max . 
‘Qualitatively, the bead is more likely to jump from one state to another 
thanks to thermal fluctuations if the waiting time is longer. Of course it also 
has fewer chances to escape from state 0 if the force pushing it toward this 
state is stronger. We did some measurements keeping the product r x / max 
constant. The results are shown in figure 7 (blue points). 

The proportion of success is clearly not constant when the product r x / max 
is kept constant, but, as expected, the higher the force, the higher P s . One 
must also note that the experimental procedure never reaches a Ps higher 
than ~ 95%. This effect is due to the last part of the procedure: since the 
force is stopped when the barrier is low, the bead can always escape from 
state 0 during the time needed to rise the barrier. This problem can be 
overcome with a higher barrier or a faster rising time Thigh- It was tested 
numerically by Raoul Dillcnscheider and Eric Lutz, using a protocol adapted 
from [6] to be close to our experimental procedure. They showed that for 
a high barrier of 8k B T the proportion of success approaches only ~ 94%, 
whereas for a barrier of 15 k B T it reaches ~ 99 %. Experimentally we define 
a proportion of success Ps force by counting the number of procedures where 
the particle ends in state 0 when the force is stopped (before the rising of the 
barrier). It quantifies the efficiency of the pushing force, which is the relevant 
one since we have shown that the pushing force is the only contribution to the 
mean dissipated heat. Measured Ps force are shown in figure 7 (red points). 
Ps force is roughly always 5% bigger than P s and it reaches 100% of success 
for high forces. 

To reach the Landauer’s bound, the force necessary to erase information 
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max 


Figure 7: Proportion of success for different values of r and / max , keeping 
constant the product r x / max ~ 0.4 pN s (which corresponds to r x v max = 
20 pm). Ps quantihes the ratio of procedures which ends in state 0 after the 
barrier is risen. Ps f orce quantihes the ratio of procedures which ends in state 
0 before the barrier is risen. 


must be as low as possible, because it is clear that a higher force will always 
produce more heat for the same proportion of success. Moreover, the bound 
is only reachable for a quasi-static (i.e. r —» oo) erasure procedure, and the 
irreversible heat dissipation associated with a finite time procedure should 
decrease as 1/r [18]. Thus we decided to work with a chosen r and to 

manually 5 optimise the applied force. The idea was to choose the lowest 
value of / max which gives a Ps force >95%. 

The Landauer’s bound k^T In 2 is only valid for totally efficient procedures. 
Thus one should theoretically look for a Landauer’s bound corresponding to 
each experimental proportion of success (see equation 1). Unfortunately the 
function In 2 +p In (p) + (1 — p) ln(l — p) quickly decreases when p is lower than 
1. To avoid this problem, we made an approximation by computing (Q )^o 
the mean dissipated heat for the trajectories where the memory is erased 
( i.e the ones ending in state 0). We consider that (Q)^o mimics the mean 
dissipated heat for a procedure with 100 % of success. This approximation is 
reasonable as long as Ps force is close enough to 100%, because the negative 
contributions which reduce the average dissipated heat are mostly due to the 
rare trajectories going against the force (i.e ending in state 1). Of course, 

5 The term “manually” refers to the fact that the optimisation was only empirical and 
that we did not computed the theoretical best / max for a given value of r. 
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at the limit where the force is equal to zero, one should find Ps force = 50% 
and (Q)- 5.0 = 0 which is different from ^^ 1112 . The mean dissipated heat 
for several procedures are shown in figure 8. 
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Figure 8: Mean dissipated heat for several procedures, with fixed r and 
different values of / max . The red points have a force too high, and a Ps force > 
99%. The blue points have a force too low and 91% < Ps force <95% 
(except the last point which has Ps force ~ 80%). The black points are 
considered to be optimised and have 95% < Ps f orce < 99%. The errors 
bars are ±0.15 k B T estimated from the reproductibility of measurement with 
same parameters. The fit (Q )~>0 = In 2 + B/t is done only by considering 
the optimised procedures. 


The mean dissipated heat decreases with the duration of the erasure pro¬ 
cedure r and approaches the Landauer’s bound k^T \\\2 for long times. Of 
course, if we compute the average on all trajectories (and not only on the 
ones ending in state 0) the values of the mean dissipated heat are smaller, but 
remain greater than the generalised Landauer’s bound for the corresponding 
proportion of success p: 

(Q)-> 0 > (Q) > k B T [In 2 + pln(p) + (1 -p)ln(l - p)] (15) 

For example, the last point (r = 40 s) has a proportion of success Ps fo rce ~ 
80%, which corresponds to a Landauer’s bound of only ~ 0.19 k B T, and 
we measure (Q )^0 = 0.59 k B T greater than ( Q) = 0.26 k B T. The manually 
optimised procedures also seem to verify a decreasing of (Q)^ 0 proportional 
to 1 /t. A numerical least square fit (Q )^0 = In 2 + B/t is plotted in figure 8 
and gives a value of B = 8.15 k B T s. If we do a fit with two free parameters 
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{Q)^o = A + B/t, we find A = 0.72 k^T which is close to k^T In 2 ~ 
0.693 fc B r. 

One can also look at the distribution of Qdrag-s-o- Histograms for procedures 
going from 1 to 0 and from 0 to 0 are shown in figure 9. The statistics are not 
sufficient to conclude on the exact shape of the distribution, but as expected, 
there is more heat dissipated when the particle has to jump from state 1 to 
state 0 than when it stays in state 0. It is also noticeable that a fraction 
of the trajectories always dissipate less heat than the Landauer’s bound, 
and that some of them even have a negative dissipated heat. We are able to 
approach the Landauer’s bound in average thanks to those trajectories where 
the thermal fluctuations help us to erase the information without dissipating 
heat. 




(a) 1 —>• 0 (b) 0 —>■ 0 

Figure 9: Histograms of the dissipated heat Qdrag-sn- (a) For one procedure 
going from 1 to 0 (r = 10 s and / max = 3.8 x 10 -14 N). (b) For one procedure 
going from 0 to 0 (r = 5 s and / max = 3.8 x 10~ 14 N). The black vertical lines 
indicate Q = 0 and the green ones indicate the Landauer’s bound k^T In 2. 
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4 Integrated Fluctuation Theorem applied on 
information erasure procedure 

We have shown that the mean dissipated heat for an information erasure 
procedure applied on a 1 bit memory system approaches the Landauer’s 
bound for quasi-static transformations. One may wonder if we can have 
direct access to the variation of free-energy between the initial and the final 
state of the system, which is directly linked to the variation of the system’s 
entropy. 

4.1 Computing the stochastic work 

To answer to this question it seems natural to use the Integrated Fluctuation 
Theorem called the Jarzinsky equality [7] which allows one to compute the 
free energy difference between two states of a system, in contact with a heat 
bath at temperature T. When such a system is driven from an equilibrium 
state A to a state B through any continuous procedure, the Jarzynski equality 
links the stochastic work W st received by the system during the procedure 
to the free energy difference A F = Fb — Fa between the two states: 

( e -/3W st ) = e -0A F ( 16 ) 

Where (.) denotes the ensemble average over all possible trajectories, and 

For a colloidal particle confined in one spatial dimension and submitted to 
a conservative potential V(x, A), where A = A(t) is a time-dependent external 
parameter, the stochastic work received by the system is defined by [16]: 

W„(x(t)} = J'^kdt' (IT) 

Here the potential is made by the double-well and the tilting drag force 6 

V(x, A) = U 0 (x,I(t)) - f(t)x (18) 

and we have two control parameters: I(t) the intensity of the laser and fit) 
the amplitude of the drag force. 

Once again, we can separate two contributions: one coming from the low¬ 
ering and rising of the barrier, and one coming from the applied external drag 

6 Since we are in one dimension, any external force can be written as the gradient of a 
global potential. 
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force. We again consider that the lowering and rising of the barrier should 
not modify the free-energy of the system, and that the main contribution is 
due to the drag force. Thus: 


W Bb = 


f*^low+ r 


' ^Tow 


-fx d t! 


(19) 


Noting that f(t = T low ) = 0 = f(t = T [ow + r), it follows from an integration 
by parts that the stochastic work is equal to the heat dissipated by the drag 
force: 

r'l iow+ r /‘^iow+ r 


W st = 


-fxdt' = 


f x d t’ = Q d 


'2iow 


rag 


( 20 ) 


The two integrals have been calculated experimentally for all the trajectories 
of all the procedures tested and it was verified that the difference between 
the two quantity is completely negligible. In the following parts, we write 
W st for theoretic calculations, and Qdrag when we apply the calculations to 
our experimental data. 


4.2 Interpreting the free-energy difference 

Since the memory erasure procedure is made in a cyclic way (which implies 
A U = 0) and AS = — k-Q\\\2 it is natural to await A F = kp,T In2. But the 
A F that appears in the Jarsynski equality is the difference between the free 
energy of the system in the initial state (which is at equilibrium) and the 
equilibrium state corresponding to the final value of the control parameter: 

ATj arzynski = -F(A(ffi na i)) — F(A(£ in j t i a l)) (21) 

Because the height of the barrier is always finite there is no change in the 
equilibrium free energy of the system between the beginning and the end of 
our procedure. Then AFj arzynsk i = 0, and we await (e~ /3H/st ) = 1, which is 
not very interesting. 

Nevertheless it has been shown [22] that, when there is a difference between 
the actual state of the system (described by the phase-space density p t ) 
and the equilibrium state (described by p® q ), the Jarzynski equality can be 
modified: 

/„-/3w st (m _ Pf q (M(X>) /3AFj arzynski (t) (00} 

p ( x ,t ) e [ZZ) 

Where (.)/ ^ is the mean on all the trajectories that pass through x at t. 

In our procedure, selecting the trajectories where the information is ac¬ 
tually erased is equivalent to fix the position x to the chosen final well 


16 



(state 0 corresponds to x < 0) at the time t = Ti ow + r. It follows that 
p{x < 0, T low + r) is directly Ps force, the proportion of success of the proce¬ 
dure, and p eq (x < 0, A(T low +r)) = 1/2 since both wells have same probability 
at equilibrium'. Then: 


2 -/3W s t(T low +r)\ 1/2 


/ ->o 


P 


S force 


(23) 


Similarly for the trajectories that end the procedure in the wrong well (state 
1 ) we have: 


2 -/3M4t(T low +r)^ _ 


1/2 


1 -P 


S force 


(24) 


Taking into account the Jensen’s inequality, i.e. (e x ) > e ( x \ we find that 
equations 23 and 24 imply: 


(Ihst)^o ^ k ® T [ ln ( 2 ) + In(Tsforce)] (e){ -s 

{w*)^ > k B T [ln(2) + ln(l - P 5force )] 1 j 

Since that the mean stochastic work dissipated to realise the procedure is 
simply: 

(khst) = Ps force (Wst)-^ + (1 — ^Sforce) (khst)^! (26) 

it follows: 


(khst) > k B T [ln(2) + Pg force hi ( Ps force) + (1 — -^Sforce) ln(l — P^force)] (27) 

which is the generalization of the Landauer’s bound for P s f orce < 100%. 
Hence, the Jarzynski equality applied to the information erasure procedure 
allows one to find the complete Landauer’s bound for the stochastic work 
received by the system. 

Finally we experimentally compute AP e fj which is the logarithm of the 
exponential average of the dissipated heat for trajectories ending in state 0: 

AP eff = - In ((e _/3Qdrag >^ 0 ) . (28) 

Data are shown in figure 10. The error bars are estimated by computing 
the average on the data set with 10% of the points randomly excluded, and 
taking the maximal difference in the values observed by repeating this oper¬ 
ation 1000 times. Except for the first points * 8 (r = 5 s), the values are very 

' A more detailed demonstration is given in Appendix. 

8 We believe that the discrepancy can be explained by the fact that the values of Qdrag ->o 
are bigger and that it is more difficult to estimate correctly the exponential average in this 
case. 
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Figure 10: Mean dissipated heat (*) and effective free energy difference (x) 
for several procedures, with fixed r and different values of / max . The red 
points have a force too high, and a Ps force > 99%. The blue points have a 
force too low and 91% < Psiorce <95% (except the last point which has 
Ps force ~ 80%). The black points are considered to be optimised and have 
95% < Ps force < 99%. 


close to k^T In 2, which is in agreement with equation 23, since Ps force is 
close to 100%. Hence, we retrieve the Landauer’s bound for the free-energy 
difference, for any duration of the information erasure procedure. 

Note that this result is not in contradiction with the classical Jarzynski 
equality, because if we average over all the trajectories (and not only the ones 
where the information is erased), we should find: 


(e~P Wst > = P Sforce (e-^‘)^ 0 + (1 - Ps force) (e"* 1 )^ = 1. (29) 


S force 


However, the verification of this equality is hard to do experimentally since 
we have very few trajectories ending in state 1, which gives us not enough 
statistics to estimate (e~ /3Wst )_ >1 properly. 

4.3 Separating sub-procedures 

To go further, we can also look at the two sub-procedures 1 —> 0 and 0 —* 0 
separately. To simplify calculations, we make here the approximation that 

Ps force = 100%. 

We can compute the exponential average of each sub-procedure: 


M M = (e- w ‘>^ 0 and M 0 ^ 0 = (e“ Wt > 0 ^ 0 (30) 
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For each sub-procedure taken independently the classical Jarzynski equality 
does not hold because the initial conditions are not correctly tested. Indeed 
selecting trajectories by their initial condition introduces a bias in the initial 
equilibrium distribution. But it has been shown [9] that for a partition of the 
phase-space into non-overlapping subsets \j (J — 1 , ■■■, K) there is a detailed 
Jarzynski Equality : 

= h 


(e~P Wat ) 


(31) 



with: 


Pi = / p(4) dicdp and pj = / p[t a ) dxdp 


'Xj 


'Xj 


(32) 


where p{t a ) and p(t a ) are the phase-space densities of the system measured at 
the same intermediate but otherwise arbitrary point in time, in the forward 
and backward protocol, respectively. The backward protocol is simply the 
time-reverse of the forward protocol. 
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Figure 11: Exponential means computed on the sub-procedures, for several 
parameters, with fixed r and manually optimised values of / max - The error- 
bars are estimated by computing the exponential mean on the data set with 
10% of the points randomly excluded, and taking the maximal difference in 
the values observed by repeating this operation 1000 times. 


Here, we take only two subsets j = {0 —> 0,1 —* 0}, defined by the 
position where the bead starts, and we choose t a = 7j ow the starting point of 
the applied force. Then we have: 




Pi^O c -0AF pff = 

Pi-> 0 


Pp->i 1 
~l/2 2 


( 33 ) 
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where i = {0,1} and -Po-h is the probability that the system returns into 
its initial state % under the time-reversed procedure (which always starts in 
state 0). Finally: 


Mi—s.0 — P 0 ^i an d — P 0 ^o (34) 

Experimental data are shown in hgure 11. 0 is an increasing function 

of r whereas M 0 ^ 0 is decreasing with r. Their sum is always close to 1 
(which is the same result that gives AP e g « k-gT In 2), and is always 

smaller than These observations are intuitive because the work is 

higher when the bead jumps from state 1 to 0 than when it stays in state 
0, and the work is higher when r is smaller. The interpretation of eq. 34 
gives a little more information. It is indeed reasonable to think that for 
time-reversed procedures the probability of returning to state 1 is small for 
fast procedures and increases when r is bigger, whereas the probability of 
returning to state 0 increase when r is smaller 9 . 

To be more quantitative one has to measure Po->i and Po-> 0 ; but the time- 
reversed procedure cannot be realised experimentally, because it starts with 
a very fast rising of the force, which cannot be reached in our experiment. 
Thus, we performed numerical simulations, where it is possible to realise the 
corresponding time-reversed procedure and to compute Po->i and Po->o- We 
simply integrate eq. 3 with Euler’s method, for different set of parameters 
as close as possible to the experimental ones. The Gaussian white noise 
is generated by the “randn” function from Matlab® (normally distributed 
pseudorandom numbers). For each set of parameters we repeat the numerical 
procedure a few thousand of times. Some results are shown in table 1 (values 
are estimated with errorbar ±0.02): 


r (s) 

/max (fN) 

Mi^o 

Pq->i 

d7 0 _).o 

^o-o 

Ps (%) 

Ps force (%) 

5 

37.7 

0.17 

0.16 

0.86 

0.84 

97.3 

99.8 

10 

28.3 

0.29 

0.28 

0.74 

0.72 

96.6 

99.3 

20 

18.9 

0.42 

0.41 

0.63 

0.59 

94 

97.1 

30 

18.9 

0.45 

0.43 

0.59 

0.57 

94.4 

97.7 


Table 1: Results for simple numerical simulations of the experimental proce¬ 
dure. 

9 Of course Po->i + -Po->o = 1- 
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All the qualitative behaviours observed in the experimental data are re¬ 
trieved, and the agreement between Mj_> 0 and P 0 -h is correct. It was also 
verified that for proportions of success < 100%, if one takes all the tra¬ 
jectories, and not only the ones where the bead ends in the state 0, the 
classical Jarzynski equality is verified: (e - ^ 1 ^) = 1. This result means that 
the small fraction of trajectories where the bead ends the erasure procedure 
where it shouldn’t, which represent sometimes less than 1 % of all trajecto¬ 
ries, is enough to retrieve the fact that AF Jarzynski = 0. 
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5 Conclusion 


In conclusion, we have realised an experimental information erasure proce¬ 
dure with a 1-bit memory system, made of a micro-particle trapped in a 
double well potential with optical tweezers. The procedure uses an external 
drag force to reset the memory of the system in one state, i.e. erase the 
knowledge of previous state and lose information. We measured the propor¬ 
tion of success of erasure procedures with different duration r and amplitude 
of the force / max . These data were used to manually optimised the procedure, 
i.e. for different r we found the lowest force which gives a good erasure of 
information. By varying the duration of the information erasure procedure 
r, we were able to approach the Landauer’s bound k-^T In 2 for the mean 
dissipated heat by the system (Q). We have also shown that (Q) seems 
to decrease as 1/r, which is in agreement with the theoretical prediction 
for an optimal information erasure procedure [1], and was later confirmed 
experimentally with a more controlled experimental system [8]. 

We have computed the stochastic work received by the system during the 
procedure, which is in our particular case equal to the heat dissipated by 
the action of the external force. We used a modified version of the Jarzyn- 
ski equality [22] for systems ending in a non- equilibrium state to retrieve 
the generalised Landauer’s bound for any proportion of success on the mean 
stochastic work received by the system. This relation has been tested exper¬ 
imentally, and we have shown that the exponential average of the stochastic 
work, computed only on the trajectories where the information is actually 
erased, reaches the Landauer’s bound for any duration of the procedure. 

We also used a detailed version of the Jarzynski equality [9] to consider 
each sub-procedure where the information is erased (1 —y 0 and 0 —» 0) 
independently. This relation allowed us to link the exponential average of 
stochastic work, computed only on a subset of the trajectories (corresponding 
to one of the sub-procedures), to the probability that the system returns to 
its initial state under a time- reversed procedure. We have shown that the 
experimental data are qualitatively in agreement with this interpretation. 
Finally, we used some very simple numerical simulations of our experimental 
procedure to compare quantitatively the partial exponential averages to the 
probabilities that the system returns in its initial state under time-reversed 
procedures. 
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Appendix 


Equation 23 is obtained directly if the system is considered as a two state 
system, but it also holds if we consider a bead that can take any position in 
a continuous ID double potential along the x-axis. We place the reference 
x = 0 at the center of the double potential. 

Equation 22 states: 

/ p -/3W st (t)\ _ P eq ( x , Ht)) r -0AF, ar ^^(t) /qep 

' h*,t) p ( x ,t) 

where is the mean on all the trajectories that pass through x at t. 

We choose t = Ti ow + r the ending time of the procedure, and we will not 
anymore write the explicit dependence upon t since it’s always the same 
chosen time. We recall that AFj arzyns ki = 0 at t = T\ ow + r for our procedure. 
We define the proportion of success, which is the probability that the bead 
ends its trajectory in the left half-space x < 0: 

Psforce = p(x < 0) = f d X p{x) (36) 

J —oo 

The conditional mean is given by: 

( e-/3Wst >* = J dW Bb p(W Bt \x)e-f* w « (37) 


where p{W st \x) is the conditional density of probability of having the value 
Wst for the work, knowing that the trajectory goes through x at the chosen 
time r. 

We recall from probability properties that: 


p(W st \x) 


p(W st ,x) 

p{x) 


(38) 


where p(W s t, x) is the joint density of probability of the value W st of the work 
and the position x through which the trajectory goes at the chosen time r. 
We also recall: 


p(W st \x < 0) 


J^dxpjW^x) 

f-oo dx P(x) 


f-oodxpjW^x) 

force 


(39) 


Then by multiplying equation 35 by p(x) and integrating over the left half¬ 
space x < 0 we have: 



dx p eq (x) 


(40) 
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Since the double potential is symmetric J 0 ^ dxp eq (x) = 
By applying definition 37 and equality 38, it follows: 

J dx J dW st p(W st ,x)e~P Wat = ^ 

Then using equality 39: 

Ps force J dW stP {W st \x < 0)e~^ = l 2 

Finally we obtain: 

/ e -pw Bt \ = V 2 

\ /a;(T low +r)<0 P s {orce 

which is equation 23 of the main text. 


( 41 ) 


(42) 


(43) 
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